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SUMMARY 


An analysis has been made of the stability of a system composed of 
a subsonic canard airframe and a canted-axis gyroscope automatic pilot. 
The analysis was based on calculations of the airframe frequency response 
and measurements of the autopilot frequency response. These responses 
were combined in accordance with the principles of servomechanism theory 
so that the stability could be determined on the basis of the Nyquist 
criterion. 


INTRODUCTION 


The Langley Aeronautical Laboratory of the National Advisory 
Committee for Aeronautics has conducted an investigation of the sta- 
bility of a system composed of a BUbsonic canard airframe and canted- 
axis gyroscope control equipment. The object of the investigation was 
to develop a method of analysis that could be used to determine the 
stability and controllability of the system and to show the application 
of this method by means of an illustrative example. 

The analysis was based on measured autopilot frequency responses 
obtained from bench tests of the control equipment and calculations of 
the airframe response. These responses were combined in accordance with 
the principles of servomechanism theory so that the system stability 
could be determined on the basis of the Nyquist criterion. The tran- 
sient response of the airframe and autopilot system to a command signal 
was determined by the use of a Fourier series technique. 
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DEFINITIONS AND SYMBOLS 


Vector frequency response: A vector which is a frequency-dependent 

function. The angle of the vector is defined by the phase relation- 
ship between the input or forcing function and the output or response , 
function and the magnitude is the ratio of the peak amplitude of the 
output function to the peak, amplitude of the input function when the 
time variation of the input is sinusoidal. 

Closed— loop response: The vector frequency response of a system 

which is sensitive to error signal. The error is the difference between 
the input signal and the output signal. 

The closed— loop response is equal to the ratio of the output signal 
to the input signal. 

Open— loop response: The vector frequency response of a closed-loop 

system with the feedback link open. The open link may be indicated by 
a broken line in block diagrams. 

The open— loop response is equal to the ratio of the output signal 
to the error signal. 

Nyquist diagram: A polar plot of the locus of the tip of the open- 

loop frequency— response vector. This diagram indicates stability of the 
closed— loop system if the critical point (—1,0° is the critical point for 
all diagrams presented herein) is not enclosed by the locus. 

Gain margin: On a Nyquist diagram, the percentage gain as measured 

by the distance from the critical point to .the point where the locus 
cuts the 180° line. 

Phase margin: The angle in degrees between the l8o° line and the 

line from the origin to the point where the locus cuts the unit circle 
on a Nyquist diagram. 

Gyroscope cant angle: Angle between the yaw axis and the spin axis 

when the spin axis lies in the plane of symmetry of the airframe. 

Follow-up ratio: The ratio of gyroscope— error angle to control- 

surface deflection. 

Z • stability axis which lies in the plane of symmetry and 

perpendicular to the relative wind and passes through 
the center of gravity 

X stability axis which lies in the plane of symmetry an d 

perpendicular to the Z-axis and passes through the 
center of gravity 



NACA TN 2004 


3 


T 

ix 

X T 

T Z 

k X 

k z 

m 

S 

c 

b 

Y 

P 

q. 

r 

3 

a 

n 

7 

e 

0 


stability axis which passes through the center of gravity 
and is perpendicular to the plane of symmetxy 

moment of inertia about the X-axis , slug-feet square 

moment of inertia about the Y-axis, slug-feet square 

moment of inertia about the Z-axis, slug-feet square 

radius of gyration about the X-axis, feet 
radius of 'gyration about the 2-axis, feet . 

mass, slugs 
wing area, square feet 
wing chord, feet 
wing span, feet 
velocity, feet per second 

rolling angular velocity, degrees per second 

dynamic pressure, pounds per square foot (H or 

pitching angular velocity, degrees per second 

yawing angular velocity, degrees per second 

angle of sideslip, positive when the relative wind is 
to the right of the longitudinal body axis, degrees 

angle of attack, degrees 

angle between the longitudinal flight path and the longi- 
tudinal principal axis of inertia, positive when 
principal axis is above the flight path, degrees (fig. l) 

climb or glide angle measured from horizontal, positive in 
climb, degrees (fig. l) 

angle between longitudinal principal axis and body axis, 
degrees (fig. l) 

pitch angle, nose-up is positive, degrees (fig. l) 
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roll angle, positive in roll to right, degrees 

yaw angle, positive when nose is to right, degrees 

elevator angle, positive when trailing edge is down, 
degrees 

aileron angle, positive when right aileron is down, 
degrees 


lift coefficient 


Off) 


variation of lift coefficient with 


angle of attack 



/d C T \ 

variation of lift coefficient with elevator angle — — 

Vd bj 


variation of pitching-moment coefficient with angle of 
attack 


© 


variation of pitching-moment coefficient with elevator 
angle 


(QA 

\Ss© / 


variation of pitching-moment coefficient with pitch ing - 

/SCmV 

angular— velocity factor / — — 

ldS2. 

\ 2YJ 

variation of rolling-moment coefficient with aileron 
angle 



variation of rolling-moment coefficient with angle of 


sideslip 
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variation of side— force coefficient with angle of 
sideslip 

variation of yawing-moment coefficient with angle of 
sideslip 

variation of yawing-moment coefficient with rolling- 



angular— velocity factor 


2VJ 


Clr 


- y n r 


Cn St 




variation of rolling-moment coefficient with rolling- 


angular— velocity factor 


r Sc 1 

CFb 

id 2 v; 


variation of yawing-moment coefficient with yawing- 


angulaj>-velocity factor 


^rb 
‘ 2V> 


variation of yawing-mioment coefficient with aileron 
angle 

variation of side— force coefficient with aileron 



angle 


/ SCy\ 

WJ 


variation of rolling-moment coefficient with yawing— 

r ac. 

angular-velocity factor ( 

h — — i 
2V/ 


D 


differential operator 


/_d_\ 
\dt/ ' 



6 


MCA TN 2004 


a) frequency, radians per second 

i = fa 


e 


DB 




phase angle, positive when output leads input and negative 
when output lags input, degrees; or error signal 


decibel gain where one decibel equals 20 times the 

logarithm to the base 10 of the absolute value of the 
ratio of the output function to the input function 


autopilot vector frequency response [%/d, 5 a /0, or 



Aq, Aj, A^ airframe vector frequency response, mode of motion 
^ indicated by subscript 

Subscripts: 

i input or forcing function, corresponding to a command 

calling for a change in heading or to a sinusoidal 
input variation 


o output or response function, for example, the system 

response to a command signal or to a sinusoidal input 
variation 


APPARATUS AND TECHNIQUE 


Airframe .— The airframe is a canard type. Thrust is derived from 
a pulse— jet— type motor which is mounted above the fuselage. A photo- 
graph of the airframe is included as figure 2. Longitudinal control is 
attained by the use of horizontal hinged elevator surfaces placed near 
the nose of the fuselage. Lateral control (both roll and yaw) is 
obtained from conventional ailerons at the wing trailing edge. Direc- 
tional stability is obtained from fixed vertical fins at the rear of the 
fuselage. The airframe has no movable control surface corresponding to 
a conventional rudder. 

Automatic pilot .— The automatic pilot is an elementary, inexpensive, 
canted-axis, free gyroscope with two sets of electrical contact pick- 
offs, one set (sensitive to pitch angle) transmitting gyroscope intelli- 
gence to the elevators and the other set (sensitive to yaw and roll angle) 
carrying signals to the ailerons. Each set of pick-offs has a dead spot 
through which the gyroscope contact must move before an error signal may 
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be transmitted to the servomotor actuating the corresponding control 
surfaces. When the dead spot is traversed in response to a disturbance, 
a signal is carried through the contacts to the servomotor, which moves 
the control surface and also actuates the autopilot pulley-and-cable 
follow-up system to center the gyroscope pick— off in the dead spot. 

The autopilot thus provides on— off control of the constant— speed servo- 
motors and, within the limitations imposed by the dead spot and the 
limitations imposed by the dynamic response of the constant— speed 
motors, causes correcting control— surface deflections proportional to 
the angle to which the airframe is disturbed. 

Remote control is provided by precession coils in the autopilot, 
which, when energized due to a remote-control sig na l, cause the gyro- 
scope to precess and thus change the space orientation of the airframe. 

A photograph of the autopilot installed in the airframe is shown as 
figure 3* 

Test apparatus .- The ratios of gyroscope-error angle to control- 
surf ace-deflection angle (follow-up ratio) were determined by mounting 
protractors on the autopilot and control surfaces of the airframe. 

An autopilot and two servomotors were mounted so that the follow- 
up ratios duplicated the ratios used in the airframe. This autopilot 
servomotor system was mounted on a table which was oscillated sinus- 
oidally in a range of varying amplitudes and frequencies. Table 
motion end servo motion were measured by NACA control-position recorders. 
These motions, as well as 0. 1-second time sig n als, were recorded by a 
galvanometer mounted on a standard NACA instrument base. The autopilot 
test setup is shown in figure 4. 

Preliminary test measurements .— Follow-up ratios were determined 
by measuring the control— surface deflections resulting from measured 
angles of gyroscope displacement in the complete airframe. In addition, 
sufficient measurements were made to determine the ratios of servomotor 
deflections to corresponding control-surface deflections. These ratios 
were used to convert servo motion (as measured in the frequency— response 
tests described subsequently) into control-surface motions. Table I 
summarizes the autopilot characteristics. 

Frequency— response tests .— The vector frequency response of the 
autopilot servomotor control system was determined by subjecting the 
autopilot to sinusoidal oscillations. The autopilot was oscillated in 
roll, pitch, and yaw. All the autopilot frequency— response data used 
in the analysis were obtained from the test setup which incorporated 
the same pulley and follow-up arrangement as the airframe. 

Inasmuch as preliminary studies of the system indicated a large 
stability margin for most conditions, frequency— response data were read 
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by an approximate method, thereby the time and labor involved in data 
reduction were greatly reduced. The accuracy of these data is therefore, 
estimated at ±10 percent. Although a more exact method was suggested in' 
reference 1, the method used herein was considered sufficient for this 
investigation . 

No attempt was made to si mu late the servo load represented by hinge 
moments while the autopilot frequency— response characteristics were 
being determined. Tests showed that when torques equal to the torques 
due to the hinge moments were applied to the servomotor, negligible 
effects were produced. 

Airframe frequency-response calculations .- The airframe vector— 
frequency-response calculations were made in accordance with the method 
outlined in reference 1 for pitch, roll, and yaw, under three flight 
conditions: 

(a) Level flight at 400 miles per hour 

(b) Level flight at 600 miles per hour 

(c) Climb angle of 3° at 250 miles per hour 

Constant speed was assumed for each flight condition,- thereby 
the calculations were simplified. 

The equations of motion used for these calculations may be found 
in the appendix. The values of airframe characteristics used in 
arriving at the solutions may be found in table II. 


METHOD OF ANALYSIS 


The entire analysis of the system stability is based on the deter- 
mination of the vector frequency response of the components. Since the 
airframe is to be stabilized in pitch, roll, and yaw, responses in these 
three planes are required. Proper manipulation of these responses anH 
proper operation upon them will allow the determination of system 
stability, of the type of transient response to be expected, and of 
recommendations for changes in system parameters. 


The derivation of open— loop and closed— loop response equations for 
the system follows: 

The error angle is denoted by e and is the angle between the 
gyroscope reference and the instantaneous airframe heading, the 



NACA.TN 2001+ 


9 


subscript ^ or 0 indicates an error due to i|r or 0. The method 
of the derivation found herein is carried through in general form in- 
references 2, 3* and 5» 

The airframe and autopilot lateral system may he drawn in' block 
diagram form as follows: 



This is a multiloop servomechanism and a succession of Ny quiet 
diagrams may be used to determine its stability (reference 2, p. 73)- 
For convenience, the figure may be redrawn to show only the inner loop 
since its stability must be examined first. 



following equations are solved: 

= ^p e 0 
to = A0 5 a 

Therefore, . 
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From a Nyquist plot of this response, the inner— loop stability may 
be determined. 

If this loop is stable, the usual, Nyquist criterion may be applied 
to the complete lateral system and, for this purpose, the system open- 
loop response is plotted in polar form. This response may be 

found more easily if the block diagram is redrawn as follows: 



The lateral system open- loop response is equal to the inner— loop 
closed— loop response 0 O times the airframe response, A^/A^. The 

response 0 o /€^ may be obtained from the following equations: 


t L 


= e t “ $o 


Then 

$o _ A p A 0 
1 + Ap 


The relationship between the open— loop and closed— loop frequency- 
response vectors is illustrated graphically by figure 5. 
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The system open— loop response is 


!o _ ApAfl 

H ~ 1 + Vfy A 0 


or 


t 0 _ A p A t 

^ " 1 + VV 


Since the lateral system closed-loop response may he used to find 
the yaw transient response, the following derivation is given: For the 

closed loop, ^ — t Q . Substituting this value of into the 

expression for the open— loop response gives 


to = VSr 

+i " 1 t Vt + V0 


The longitud inal system may he drawn in block diagram form, as 
follows : 



This is a single-loop servomechanism and the open-loop response 
is e 0 /e, which equals ApA Q . The closed-loop response Q 0 /e ± is equal 

to ApA ^ — . Again, this relationship is illustrated graphically by 

1 + ApA 0 

figure 5. The airframe-autopilot system variations which were considered 
in the analysis are presented in table III. 
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RESULTS AND DISCUSSION 
Longitudinal Characteristics 


Autopilot frequency response .— Experimentally determined curves 
showing the response of the autopilot to variahle— frequency sinusoidal 
oscillations of different amplitudes are presented in figure 6. These 
curves were obtained by reading the film record to an accuracy of 
only ±10 percent since a preliminary analysis indicated that the system 
would have a large margin of stability. 

/ 

Airframe frequency response .— The calculated curves showing the 
response of the airframe to sinusoidal elevator deflection of variable 
frequency are presented in figure 7. The flight conditions shown were 
included so that the system response during a typical climbing condition - 
and during high-speed flight might be determined. 

Open— loop response .— The autopilot response and airframe response 
were combined in the manner shown under method of analysis to form the 
open— loop response. 

These combined data are presented in the form of Nyquist diagrams. 

A Nyquist diagram is a polar plot of the vector frequency response. 

The curve plotted is the locus of the tip of the response vector. This 
is a compact form of presentation which allows the designer to plot all 
significant data in one curve on one chart. These diagrams indicate 
system stability if the appropriate critical point is not enclosed, 
while a measure of the damping may be obtained from the nearness of the 
locus to, the critical point. This nearness may be expressed as the gain 
margin and phase margin. A further explanation of this form of presen- 
tation, as well as a rigorous proof of the criterion for stability, may 
be found in references 2, 3, and 6. 

For all Nyquist diagrams presented herein, the scale of the radial 
coordinate is linear between values of zero and one, whereas for values 
greater than one a logarithmic (decibel) scale is used since a greater 
range of values may be plotted. 

A typical open— loop response for the longitudinal mode .of motion 
is shown in figure 8. In table IV. the gain and phase margins for the 
various conditions considered are summarized. 

Transient response .— The longitudinal transient response is deter- 
mined from the closed— loop response and the response of the closed— loop 
system to the first half— cycle of a square wave. Such a response will 
very closely approximate the response to a step function. Determination 
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of step-function response involves the response to a Fourier integral, 
a procedure -which is not only laborious but which presupposes an 
analy tical expression for the response as a function of frequency . 
However, the response to a square wave may be determined by operating 
on the vector closed— loop response (which may be obtained wholly or 
partly by experimental means) with a Fourier series expression for a ' 
square wave and summing the result, a procedure which is relatively 
simple. (See reference 7*) 

By this method, utilizing the closed-loop response (fig. 9) and- the 
response to a Fourier series, the transient response to a command 
input signal for the cruising conditions was found. (See fig. 10.) 

The response is oscillatory and underdamped. Since the phase and gain 
margin are somewhat smaller than recommended in reference 2, page 162, 
a transient response of this nature may be expected. 

If the ratio of gyroscope-error angle to control— surface-deflection 
angle is changed to a value of 7:1, the open-loop response curve shown 
in figure 11 is obtained. Phase and gain margin values of 58° and 
74 percent will be realized under these conditions. These values are 
closer to the recommended margins and a much more desirable transient < 
response, as shown in figure 12, is obtained. Table IV s ummar izes the 
gain and phase margin for various conditions considered. 

On the basis of a comparison of the transient responses shown in 
figures 10 and 12, the ratio of gyroscope-error angle to elevator- 
deflection angle should be 7:1. This modification should be made by 
changing the linkage between the servomotor and control surface while 
the ratio of gyroscope-error angle to servomotor-deflection angle is 
held to its present value. This change will insure the retention of the 
present autopilot frequency— response characteristics. 

Lateral Characteristics 

The autopilot vector-frequency-^ response curves for yaw and roll are 
presented in figure 13.' The autopilot responses for roll were determined 
for amplitudes of ±1.5°, ±3.94°, and ±5.43°. For yaw the autopilot 
responses were determined for amplitudes of il.3°, i4.07°, and —5 • 92 ° • 

The amplitudes of table oscillation in yaw differed from those in roll 
because an approximation, rather than a duplication, of the oscillating- 
table amplitude settings was expedient. 

In the discussion to be presented and the longitudinal analysis 
previously presented, the experimental amplitude values were rounded off 
to 1°, 4°, and 6°. The assumption was made that the 1°, 4°, and 6° fre- 
quency responses differ from the experimental values by a negligible 
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amount. This assumption is justifiable since the accuracy of the 
experimental autopilot frequency— response curves is ±10 percent and 
since, in all cases, only sma ll changes in the frequency— response curves 
occur for the differences "between the experimental amplitudes and the 
rounded figures. 


The airframe lateral curves for a number of conditions are shown 
in figure 14. These curves represent actual physical variations which 
might have occurred or which might have produced advantageous or dis- 
advantageous results in the over-all system stabilization. For instance, 
the angle of the principal axis of inertia migjht have been in error by 

the amount considered, or, if it had seemed advisable, (W might have 

°a 

been changed by modifying the aileron trim position. These variations 
did not have an appreciable effect on the over-all system stabilization; 
therefore, the normal values were used for the analysis. 


The phase angles of the airframe lateral frequency response were 
modified by subtracting l80° from the calculated values, since the sign 
convention employed in the calculations resulted in phase angles l80° 
out of phase with the true physical phase angle. The airframe lateral 
vector frequency responses presented herein are uncorrected. The 
corrected values were used in the analysis and for all lateral Nyquist 
plots. 

In order to conduct an analysis of the lateral control system, the 
system must be considered as a multiple— loop servomechanism since the 
control system is responsive to both yaw and roll and since the airframe 
must roll in order to yaw. In order to use the simple Nyquist test, the 
stability of the inner loop must first be examined. When this loop is 
proved stable, the simple Nyquist criterion may be applied to the complete 
system. A further explanation may be found in reference 2, page 73 • 


Both the yaw and roll vector— frequency— response characteristics 
were used to determine stability. The closed— loop response for yaw was 
used to determine the yaw— system transient response by the method out- 
lined for the longitudinal system. 

Inner— loop stability .— A typical Nyquist diagram of the lateral 
system inner loop is shown in figure 15. In addition, a summary of the 
gain and phase margins for this loop for the various flight conditions 
may be found in table IV. An ample margin of stability is available for 
nl 1 cases except level flight at 600 miles per hour. In 600-mile— per— 
hour level flight at 1° amplitude, the inner loop is actually unstable. 
For 6° amplitude and 600-mile— per-hour level flight the margins are 
rather small . 
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System stability.— In the case where the inner loop of a mnltiloop 
servo system has been shown to be stable * the simple Nyquist test con- 
cerned with the encirclement or intersection of the critical point by 
the system open— loop response may be used to determine system stability . 
This criterion has been applied to the 250-mile— per-hour climb condition 
and to the 400-mile-per-hour level-flight condition, and ample stability, 
as indicated by large gain and phase margins, is assured in both cases. 
These margins are s ummar ized in table XV and a typical system open- 
loop response curve is presented as figure l6. 

In the case where the inner loop of a multiloop servo system is 
unstable, system instability is not necessarily indicated. However, by 
the use of a succession of Nyquist diagrams the system may be found to 
be stable or unstable. A more complete explanation of this case may be 
found in reference 2, page 735 reference 6, page 158; or reference 3, 
page 172. This general Nyquist criterion has been applied to the 
600-mile— per-hour level-flight condition and the lateral system has been 
found unstable for the 1° oscillation amplitude. In addition, the inner 
loop has insufficient phase and gain margins for the 6° amplitude 
condition. 

Therefore, a steady oscillation with an amplitude between 1° and 6° 
occurs in roll and a rather underdamped rolling motion results when a 
command for yaw heading change is given to the system. This instability 
qnfl underdamped rolling motion constitute undesirable system operation 
and, for this reason, 600-mile— per-hour level flight would not be 
recommended. 

Yaw transient response .- By use of the closed-loop response, as 
shown in figure 17, and a Fourier series expression for a square wave, 
the transient response (to a command signal input calling for a yaw. 
heading change) was found. This response is shown in figure 18 for 
the 400-mile— per-hour level— flight condition.. 

Since the response of the system in yaw was slow, a somewhat more 
detailed examination of the factors influencing this response is made. 

The autopilot installation functions in a manner to produce zero aileron 
deflection when the angle of bank is equal to the instantaneous yaw 
heading change. This function means that in flight the aircraft will 
maintain an angle of bank equal to the instantaneous yaw heading change. 
It was known that the proportionality between the angle of bank and the 
instantaneous yaw heading change could be influenced by the cant angle 
of the gyroscope. 
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In order to relate this fact to the transient response, the tran- 
sient response was approximated by a step-by-step method in which the 
bank angle was held equal to the instantaneous yaw heading change in 
the following expression: 


= VAty 
57. 3g tan 0 


where At is time increment in seconds and g is acceleration due to 
gravity. Figure 19, which shows the use of the preceding expression 
for a level— flight condition of 400 miles per hour, was obtained by 
assuming a heading of 6°. Therefore, the initial bank angle was 6°. 
Substituting 0-6° and Aty = 1° and solving for At ' results in 
the first point of the. plot. Initial, conditions for the next step, 
since the desired heading change is then 9 °, are 0=5° and At = 1°. 
The equation is then solved for At and the next step carried out in 
a similar manner. The time at any point is the sum of the preceding 
time increments. The basic approximation assumed in this method is 
that the turn is caused by the horizontal component of the lift vector 
due to the angle of bank. 

The curve shown in figure 19 approximates the transient response 
obtained by the Fourier series method shown in figure 18. Also shown 
in this figure is a step— by— step approximation of the transient obtained 
when the bank angle is held equal to three times the instantaneous yaw 
heading change. The time required to reach the new yaw heading is 
about one— third of the time required for the present system. It seems 
then, that if the yaw transient response must be modified/ such a change 
could be made by adjusting the gyroscope cant angle. The system sta- 
bility could then be analyzed by the method used for the present sta- 
bility determination. • 


CONCLUDING REMARKS 


From the results of an investigation of the stability of a system 
composed of a subsonic canard airframe and a canted-axis gyroscope 
automatic pilot, the following procedure appears to be a suitable 
method of analysis of automat ically stabilized and remotely controlled 
aircraft : 

1 . Determine the frequency— response characteristics of the airframe 
and control components either .by calculation or. by experimental means. 
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2. Combine the component frequency responses in the proper man ner 
for the system considered' and plot a Nyquist diagram, or series of 
Nyquist diag rams as required, to determine the system stability. 

3. Adjust the system frequency response to the desired charac- 
teristics on the basis of the Nyquist diagrams. 

4. Determine the transient response of the adjusted system to 
insure tha t the adjustments have resulted in the desired response. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., October 20, 1949 
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APPENDIX 

AIRFRAME EQUATIONS OF MOTION AND SAMPLE CALCULATION FOR AIRFRAME 
FREQUENCY-RESPONSE CHARACTERISTICS IN PITCH 


The frequency— response characteristics in pitch of the airframe may 
he calculated ty solution of the following' differential equations of 
motion if the forward-speed effects are neglected: 


^5T-3qSc 


D 2 


- Cm^ 


2V 






/ mVD _ °L tan 7 V + 
\57.3qS 57.3 / 


mVD 

57.31S 


V (i) 


— 


/ C T tan y\ 


K” 57.3 )J 

* = C % e S e. 




where 



The derivative due to the lag in dowrjwash at the wing due to the 
tail was omitted from the pitching-moment equation because little infor- 
mation is available for downwash effects on a canard airframe. All sta- 
bility derivatives are in degree measure. For the level— flight con- 
dition, the solution of these equations for 0/5 e may be obtained by 

setting up the determinants 


JL 

&e 


Ct (- . ml/T) - - Cx \ 

^e \ 57-3qS L a/ 


IyD 2 

57-3qSc 





mVD 


mVD 

57.3qS 



( 2 ) 


57.3qS 



.NACA TN 2004 


19 


Numerical values, taken from table II, are substituted in the deter- 
minants at this point to facilitate the solution. When these values are 
substituted, the determinants become 


0.0215 * 0.014 

/ 0.004 (— 0.0809D - 0.0825) 

(0.00119D 2 + O.OOO383D) 0.014 

0.0809D , (-0.0809D - 0.0825) 


(3) 


By making the substitution D - ico and simplifying equation (3), then 


_9 — 0 . 00174ia> - 0.00184 

‘ 8 e 0 . 000096 i(u^ + 0.000129co 2 - 0.00084 ioo 


This expression may be written in the vector form: 


— = A + iB 
Se 


where 


E = \/a 2 + B 2 


and 



The quantities E and e are then plotted against o> to give the 
airframe vector-frequency— response curves. 
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The lateral equations of motion for a rudderless airframe, which 
may he derived from the equations given in reference 8 , are as follows: 


[s ' T.BqS^ ^ 2003 ^ + ~ °nr ^ + 

— (k 2 ~ ky 2 1j(cos q sin t|)D 2 - C 

[57 • 3qSb V Z ml »p 2V 


0 - C 3 = C 5 

* V a 


[57.3q.Sb 


2 cos 2 t] + k_ 2 sin 2 q)D 2 — C, -k. D 
' Z 7 ip 27 


0 + 


2k— ( k _2 _ k 2)( 

[57 . 3qSb V Z X / v 


cos T] sin rj )!) 2 




= °l6 a 5a 


( 


mVD °L 


57.3qS 57-3 


tan 7 


y\i - (— — V + - - Cy V = Cy 5 

/ \57.3/ \57.3gS ^3/ ^5 a a 


These equations are solved separately for 0/8 a and ^ /S a hy the 
method of determinants. 
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TABLE I 

CHARACTERISTICS OF THE CONTROL STSTEM 



Longitudinal 

Lateral 

Pick— off dead spot, deg 

1.3 

2.3 

Total gyroscope angular range between 
limit stops, deg 

48 

42 

Static follow-up ratio (error angle to 
control— surface angle) 

3.5:1 

2.5:1 

Total control— surface angular 

range, deg . . . . . . 

16.2 

22 

Maximum rate of control— surf ace travel, 
deg per sec ... 

11.7 • 

19.4 
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TABLE II 

AIRFRAME CHARACTERISTICS USED IN CALCULATING FREQUENCE RESPONSE 


l y, slug-ft 2 

l z , slug-ft 2 

Iy, Slug-ft 2 

k x* 

k z> * 

m, slugs 

S, sq ft 

c, ft 

b, ft . . 

V, level flight at 400 mph, ft per sec 

C^, level flight at 400 mph 

a, level flight at 400 mph, deg 

q, at 3500 -ft altitude and 400 mph, lh per sq ft 

c Ia> per deg . 

Clg^, per deg 

Cnia/ P©r deg . 

, per deg 

■^mq* per deg 

7 , deg 

e, deg 



• • 52 

. 1080 

. 975 

. 0.93 

. 4.23 
. 60.5 
. . 21 
. 1.87 

. 11.5 

. 586 

, 0.255 

• 3.09 

. 36k 

. 0.083 

. o.oo4 

-o.oi4 

. 0.022 

. -0.24 

0 and 3 
. 3.08 


q, deg . . . 

Cjj- , per deg 
°a 

C 2p , per deg 

0y p , per deg 

Cn p , per deg 

1 Cn p , per deg 

1 C *7 , per deg 
P 

•^Cn r > P er de 8 
c n§ a J P 0r de 8 
CY 8a .... 
^ r , per deg 


• • p 0, and - 

- 0.001 

- 0.002 

-0.016 

0.004 

.... —0 . 00007 

.... - 0.0112 

.... - 0.0122 

0.0002 and -0.0002 

0 

0.0003 


These values were estimated, whereas other stability derivatives were 
obtained from unpublished wind-tunnel data. '\naca~ 





UTOPILOT SYSTEM VARIATIONS INVESTIGATED 
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TABLE- IV 

SUMMAHI OF GAIN AND PHASE MARGINS OBTAINED FROM THE OPEN-LOOP RESPONSE 



Autopilot 

Gain 

Phase 

Flight condition 

input 

amplitude 

(deg) 

margin 

(percent) 

margin 

(deg) 

Longitudinal results; 

follow-up ratio a 3.5:l 



1 


100 

2 6 

Climb at 250 miles per hour 


±4 

78 

2 6 


L±6 

77 

36 



r±i 

100 

20 

Level flight at 400 miles per hour 


±4 

50 

19 


[±6 

58 

25 



[±1 

100 

21 

Level flight at 600 miles per hour 


±4 

20 

10 


1*6 

30 

1'6 

Longitudinal results 

; follow-up ratio 7*1 



j 

~±r 

90 

93 

Climb at 250 miles per hour 

( 

^±6 

90 

93 



~±l 

100 

60 

Level flight at 400 miles per hour 


±h 

74 

60 


~±6 

. 8o 

68 

Level flight at 6 00 miles per hour 

< 

f*L 

1*6. 

H 

ON O 

VJ 1 o 

^5 

55 

Lateral results; 

inne 3 >-loop response 


Climb at 250 miles per hour 


r«. 

55 

70 

70 

75 

Level flight at 400 miles per hour 

i 

r*i 

l*. 

26 

40 

35 

39 

Level flight at 600 miles per hour 

j 

r±i • 

L±6- 

b -17 

9 

-30 

7 

Lateral results; system open-loop response 

Climb at 250 miles per hour 

■< 

[S 

85 

88 

85 

85 

Level flight at 400 miles per hour 


[3 

90 

89 

88 

85 

Level flight' at 600 miles per hour 


[3 

91 

88 


a CSrroacope-exTor angle FollOTMlp ratlo . 

Control— surface deflection 

^Unstable condition. 
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Figure 2.- Canard airframe with rear hatches removed. 
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installation characteristics. 
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Figure 4.- Autopilot test setup showing oscillating table and recording instruments. 





Figure 5.- Graphical relationship between open-loop and closed-loop 

oc 

frequency- response vectors. Closed-loop response = — . 
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Figure 6..- Autopilot longitudinal vector-frequency- 
response. characteristics.. 
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Figure 8.- Nyquist diagram showing the longitudinal open-loop 
frequency response for ±4° autopilot input amplitude and 
400-mile-per-hour level flight. Ratio, of gyroscope-error 
angle to elevator angle equals 3-5 to 1. 
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Figure 10.- Longitudinal transient response of the airframe and autopilot system to a command signal 
calling for a sudden change of attitude of U° during 400-mile-per-hour level flight. Ratio of 
gyroscope-error angle to- elevator angle equals 3*5 to 1. 





NACA TN 2004 


39 



-90 


Figure 11.- Nyquist diagram showing the longitudinal open- loop frequency- 
response for ±4° autopilot input amplitude and 400-mile-per-hour 
level flight. Ratio of gyroscope-error angle to elevator angle 
equals 7 to 1. 






System response 
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(a-) 1° attitude change 


(b) 4° attitude change 


NACA -r^ 


T/me y t } v5ec 

(c) 6° attitude change. 


Figure 12.- Longitudinal transient response of the system to a command 
signal calling for a sudden change in attitude during too -mile -per- 
hour level flight. Ratio of gyroscope-error angle to control- surface 
angle equals 7 to 1. 
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Frequency, uj, racf/ans per sec 


(b) Yaw response. 


Figure 13 ■- Concluded. 









level flight. 
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Figure 15.- Nyquist diagram showing stability of the lateral inner loop. 
Autopilot input amplitude = ±6°; 400-mile-per-hour level flight; 

i) = 0°; Cng = 0.0002. 
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Figure l6.- Nyquist diagram showing lateral system stability. Autopilot 
input amplitude = ±6°; 400-mile-per-hour level flight; t) = 0°; 

C = 0 • 0002 • 





frequency/, 6J; racf/cws per *sec 

Figure 17.- Yaw closed- loop vector frequency response. Autopilot 
amplitude = ±6°; lOO-mile-per-hour level flight; q = 0°; 

On* = 0.0002. 
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heading of 6°. 400-mile-per-hour level flight; t) = 0°; = 0.0002, 
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Figure 19* - Approximate yaw transient response due to a command signal (obtained by step-by-step 
calculation) and showing the effect of gyroscope cant-angle variations. 


